Effective field theory at moderate finite 

density 



Boris Krippa 

Department of Physics and Astronomy, Free University of Amsterdam, 
De Boelelaan 1081, 1081 HV Amsterdam. 



Effective field theory of the in-medium nucleon-nucleon inter- 
action is considered. The effective range parameters are found to 
be of a natural scale. The low density limit is discussed both in 
perturbative and nonperturbative situations. In the nonperturba- 
tive case the attractive character of the nucleon-nucleon interac- 
tions in the 1 S'o channel leads to the nuclear superfluidity which 
is analyzed in the framework of the renormalization group. The 
numerical values of the corresponding energy gap are in agree- 
ment with results obtained by more more traditional many-body 
techniques. The S-wave part of potential energy per particle is 
calculated for different values of nuclear density. The role of pion 
and many-body effects is discussed. Problems and challenges in 
constructing the chiral theory of nuclear matter are outlined. 
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1. Introduction. 

Recently, there has been much interest in applying the Effective Field The- 
ory (EFT) M methods to study nuclear interactions. The key element of 
EFT is the counting rules allowing for a consistent expansion scheme with 
predictable theoretical errors. Another important ingredient of any EFT is a 
separation of scales so that the high frequency modes can be integrated out and 
manifest itself indirectly in the Low- Energy- Effective Constants (LEC's). In 
addition, the symmetries of the underlying "microscopic" theory put important 
constrains on the interaction terms in the effective Lagrangian, making them 
less arbitrary as compared to the sometimes more traditional phenomenolog- 
ical description. In the context of nuclear matter such a description implies 
the use of some phenomenological potentials leading, in general, to the differ- 
ent effective in-medium interactions and bringing uncontrollable uncertainties, 
related to the off-shell ambiguities. These ambiguities can, at least in prin- 
ciple, be eliminated in the EFT approach by using the freedom to work in 
the different low-energy representations of the underlying QCD. In the other 
words EFT may be useful to get rid of the long-standing problem arising in the 
more standard treatments of hadron interactions and related to the necessity 
of utilizing of the phenomenological form-factors to regulate the divergent in- 
tegrals and parameterize the complicated internal structure of hadrons. There 
is nothing wrong to use them in the calculations based on phenomenological 
models. However, the problem is that the underlying field theory must then be 
nonlocal making the whole issue of constructing the entirely consistent field- 
theoretical treatment of the low-energy hadron interactions extremely difficult 
to deal with. Moreover, when the problem under consideration is the interac- 
tion of hadrons with the external electromagnetic fields, the straightforward 
use of form-factors may lead to the conflict with gauge invariance. All these 
features make the reanalysis of nuclear matter in the framework of EFT use- 
ful and promising. The main element of nuclear dynamics is the elementary 
nucleon-nucleon interactions. Because of the large scale involved in this process 
the EFT treatment of the NN-interactions in vacuum must be nonperturba- 
tive. In the framework of EFT the problem was first considered by Weinberg 
@ who suggested to apply chiral counting rules to the set of all possible irre- 
ducible diagrams at given order which are then to be iterated in all order by 
solving Lippmann-Schwinger equation. In somewhat different approach pro- 
posed by Kaplan, Savage and Wise (KSW) |J only certain subclass of lowest 
order diagrams was summed up to all orders whereas the rest was treated as 
a perturbation so that KSW counting is applied directly to the scattering am- 
plitude. The analysis |0| of the NN (pionless!) scattering amplitude based 
on the Wilson renormalization group has demonstrated that power counting 
indeed closely follows that which was found by KSW. 
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The extension of these developments to the case of finite density is by 
no means straightforward [||]. There are several points which must be taken 
into consideration to formulate the consistent EFT at finite baryon density. 
First, being unnaturally large in vacuum the NN S-wave scattering lengths 
get significantly reduced in nuclear matter due to Pauli principle and many- 
body interactions. Second, the way of formulating the counting rules allowing 
for the consistent treatment of the many-body effects, such as the ring and 
ladder diagrams must be found. Third, the long-range part of the NN effective 
interaction related to the pion exchange in nuclear environment should be 
taken into account. In addition, the three and four nucleon interaction terms 
must be included in the EFT Lagrangian. 

Some recent studies of of the systems with the finite fermion density in the 
context of EFT have been focused on constructing of the chiral expansion for 
the dilute interacting Fermi-gas and mean-field effective chiral Lagrangian 
0. The issues of the in-medium counting rules || and nuclear binding in 
chiral limit || were also examined recently. In this paper we focus on the 
EFT approach to construct the analog of the Brueckner G-matrix and try 
to relate the in-medium LEC's and those obtained in the free space fit. It 
seems to be quite a nontrivial problem to establish such relations since the 
LEC's entering the effective potential of the in-medium NN interactions may 
be density dependent. 

Any EFT expansion has a form in k/A, the ratio of some generic mo- 
mentum involved and a characteristic short range physics mass scale. So it 
is important to establish a relevant expansion parameter. The analysis of the 
different phenomenological models indicates that the nonrelativistic mean field 
approximations, successfully describing the bulk nuclear properties at low and 
moderate energy scale, encounter difficulties when applied to the nuclear phe- 
nomena with a characteristic mass scale around 600 MeV (electromagnetic 
form-factors at relatively large momentum transfer, for instance). This is the 
region where the effect of the short range correlations (SRC) becomes impor- 
tant. Thus, it is plausible that the typical scale where the EFT description 
of nuclear matter with pointlike nucleon-nucleon interaction can be close to 
the inverse range of the SRC. One notes that the addition of the long-ranged 
pion effects should not significantly change this estimate. As the typical nu- 
clear momentum scale is of order the Fermi-momentum pp ~ 260-270 MeV the 
relevant expansion parameter is close to 0.3 - 0.4. This admittedly crude esti- 
mate can be improved if the influence of the nuclear many-body environment 
is taken into account. 

All essential features of the nucleon-nucleon interactions can, in principle, 
be obtained from the effective chiral Lagrangian 
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We consider the simplest case of the NN scattering in the 1 So state and as- 
sume zero total 3-momentum of NN pair in the medium. The inclusion of the 
nonzero total 3-momentum does not really change anything qualitatively and 
only makes the calculations technically more involved. The G-matrix, describ- 
ing the dynamics of the in-medium interactions satisfies the BG equation 
The G-matrix is given by 



Here ei and e 2 are the single-particle energies of the bound nucleons. In 
nuclear medium bound nucleon acquires the effective mass slightly different 
from that in free space. In a fully self-consistent calculations the nuclear 
mean field and chiral G-matrix should, in principle, be treated using the same 
Lagrangian. However, it requires the inclusion both pion and many-body 
effects and establishing the power counting rules at finite density. In this paper 
we explore an easier approach and use the in-medium value of the nucleon 
mass close to the accepted one in the phenomenological nonrelativistic [|10j 
mean field treatments. We used the value M = 0.75Mo in the calculations 
at pf = 1-36 fmT 1 , where we denote Mo the nucleon mass in vacuum. In 
order to check the sensitivity of the results obtained to the different choices 
of effective nucleon mass the values of M have been varied from 0.7 M to 
0.8M . The higher order changes of the results were found for the potential 
energy per particle so we will keep the value of M fixed and discuss the results 
corresponding the choice M = 0.75M . 

This paper is organized as follows. In the next section we consider an 
exactly solvable model to estimate the in-medium effective range parameters. 
Then the generalization of the effective nucleon-nucleon Lagrangian to the fi- 
nite density case is discussed. In the section 4 we consider the low-density limit 
in the case when the perturbative treatment is possible. The issue of superfluid 
nuclear gap is discussed in the section 5. Section 6 offers a brief summary of 
results for a potential energy of nuclear matter and some concluding remarks. 

2. Exactly solvable model. 

In order to get an idea about the typical size of the in-medium effective 
range parameters which, in turn, determine the typical scale of the problem it 
seems reasonable to exploit some exactly solvable model with the parameters 
adjusted to describe the experimental data or empirical results. Moreover, the 
in-medium scattering amplitude obtained in such solvable model can further be 
used to extract the values of the effective couplings in the corresponding EFT. 



G(p',p) = V(p',p) + M 



M(e 1 i P ) + e 2 (p'))-q- 



% ~Pf) 



aG(q,p), (2) 
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In some sense this is similar to the procedure usually adopted in the vacuum 
case. Firstly, one should determine the LEC's from, say, phase shifts and then 
use these LEC's to calculate the other observables. Unfortunately, there is no 
such thing as phase shift analysis in nuclear matter. One could instead rely on 
some solvable model with the parameters adjusted to describe nuclear data to 
extract the numerical values of the effective couplings. However, there is an 
important difference between vacuum and in-medium cases. In the former the 
corresponding LEC's are extracted from the model independent experimental 
data whereas in the latter the in-medium LEC's are determined from the model 
dependent scattering amplitude which may, in general, differ for the phase- 
equivalent NN potentials. In this respect it would be desirable to relate the 
vacuum and in-medium effective couplings to be able to use the experimental 
data on the NN scattering in vacuum. We will discuss this and related issues 
in more details below. 

We use a model with a simple separable interaction 



v = -\\v)(v\ 



with the form factors 



One can easily get 



rjip) 



(p 2 + /3 2 )V2 



T(k,k) 



V(k,k)~ l 



I -Mq J 



dqq 2 V(q,q) 
An 2 k 2 - q 2 



The analytic solution for the corresponding T-matrix is straightforward 



1 



T(k,k) 



V(k,k) 



-i 



1 - M, 



dqq 2 V(q,q) 
An 2 k 2 — q 2 



(3) 



(4) 



(5) 



(6) 



The G-matrix being the solution of the BG equation for the separable inter- 
action takes the form 



G(k,k) = -rf{k) 



A" 1 + £ / rfrf 



% - PF)rf(q) 
k 2 — q 2 



(7) 



The empirical value of the potential energy per particle is ~ -16 MeV and can 
be reproduced with 

A = 1.95 /3 = 0.8 fin (8) 



These parameters being substituted to the G-matrix in the low-momentum 



limit G(k — > 0) lead to the estimates a r 



0(1), where a m and r n 
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are the in-medium analogs of scattering length and effective radius. So one 
can conclude that the in-medium effective range parameters are of a natural 
size which is roughly given by the value of the Fermi-momentum at nuclear 
saturation pp ~ 1.36/m -1 . This allows us to avoid some part of the difficul- 
ties typical for the EFT for the nucleon-nucleon interaction in vacuum such 
as necessity of a special treatment for the system with the unnaturally large 
scattering length. Of course, there are many other specific points typical for 
the finite density EFT and making its consistent yet practical realization quite 
nontrivial. Some of them will be considered in what follows. As already men- 
tioned above the in-medium phenomenological amplitude may depend on the 
model of the NN interaction used. To check whether the freedom to choose 
among phenomenologically equivalent models of the NN interactions may give 
rise to the significant changes of our estimates of the in-medium effective range 
parameters we carried out the calculations with the exponential parameteriza- 
tion of the form-factors entering the separable potential. The values of a m and 
r m found are as well of a natural size ~ 0(1) and its numerical values are close 
to those obtained with the "square root" fit. So it seems that the conclusion 
about a natural size of the the in-medium effective range parameters is quite 
robust and practically model independent. Physically it looks quite natural. 
Firstly, the integrand in the BG equation does not have a pole irrespectively 
of the explicit form of the kernel. Secondly, the strength of the shallow virtual 
nucleon-nucleon bound state gets significantly reduced in nuclear medium be- 
cause of interaction with nuclear mean field leading to the inevitable decrease 
of the scattering length. 

3. Effective Lagrangian treatment. 

The solution of the BG equation with the NN potential extracted from the 
effective Lagrangian is similar to that in the vacuum case [II] and given by 



f^W" 1 ) 2 /<*,„), (9) 



G(k, PF ) C„ + Clh W + * 2 C 2 (2 - C,Is(f F )) 
where we the loop integrals are 



Upf) = / dqq n ~ l 9(q - p F ). (10) 



(27T) 

and 



I(Kpf) = ^J dq k2 _ n2 . (11) 



These loop integrals are divergent and therefore the procedure of regulariza- 
tion and renormalization must be carried out. Note that the issue of the 



6 



nonperturbative renormalization is quite a subtle problem. In contrast to the 
standard perturbative case where the usual field theoretical methods can be 
used to regularize the given divergent graphs and then renormalize the bare 
coupling constants, in the nonperturbative situation the renormalization of 
the whole integral equation must be carried out. In the case when the an- 
alytic solution for the scattering amplitude can be obtained as, for example, 
in the pionless nucleon-nucleon EFT, the renormalization of the amplitude is 
a rather straightforward procedure. However, if the explicit solution is not 
possible (this is the case for the realistic NN forces) then the special care is 
needed to perform the renormalization in a consistent way In this paper 



we follow the procedure used in Ref. |13[ to renormalize the effective NN am- 



plitude in vacuum. We subtract the divergent integrals at some kinematical 
point p 2 = —fi 2 . After subtraction the renormalized G- matrix takes the form 

1 1 M p F - k p F -ifi 

+ iri kln -r-^ ln —^U ( 12 ) 



G(k,p F ) C (/j>,Pf) + c 2k 2 C 2 {[i-,PF) 47r p F + k ' p F + i^ 

where the couplings Co(2) should now be interpreted as the renormalized quan- 
tities depending on some renormalization scale /i. It is easy to see that in the 
Pf — > limit the vacuum chiral NN amplitude is recovered. The \x depen- 
dence of LEC's is governed by the renormalization group (RG) equations. We 
require G-matrix to be a subtraction point independent quantity. Applying 
d/dfi to the expression for the G-matrix and setting dG/d/i = one can get 
the following RG equations 

SC ^- Pf) = CoC l M ( I 1 "" , + arctan(-)). (14) 
dfj, rr 2 pp + n 2 p F 

In the limit pp — > these equations transform to the ones derived by Kaplan 
et al. M. The solutions of these RG equations are 



Co (a* Pf) ~ Cq(/j,q,pf) ,^gs 

1 + i^(^ arctan (f^) ~ ^arctan(-^))C {^ Q) pF) 

and 



C 2 (n,p F ) = C 2 (^ , Pf) 



C (/j>o,Pf) 



(16) 



In order to determine the scale dependence of the effective couplings one needs 
to define its boundary values at some kinematical points. To extract these val- 
ues we equate the EFT and phenomenological expressions for the G-matrix 
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at p = Pf/2 and p = Pf/3. This is somewhat similar to the procedure used 
to get the values of the vacuum LEC's when the fit is usually done using 
the experimental phase shifts within some kinematical region. One notes 
that when pp — > the above expressions for the effective couplings reduce 
to these obtained in Ref. || The assumed value of the Fermi-momentum is 
Pf = 1.36/ra -1 . If the value \i = is chosen as a subtracting point we find 
C*o(/i = 0,pf) = — 1.88/m 2 in LO. For the lowest order coupling fixed at NLO 
we get C (/i = 0,pf) = — 2.67/m 2 and C 2 (/i = 0,pp) = 0.85/m 4 . Substitut- 
ing these values in the RG equations we get C (/i = m w ,pp) = 2.35 fm 2 and 
C 2 (/i = m n ,pp) = 0.64/m 4 . The obtained effective coupling C at fi = is 
fairly close to its vacuum value ||. Thus one can conclude that the density 
dependence of the in-medium LO couplings is rather moderate, provided that 
fi ~ p F . It means that at densities smaller than the normal nuclear one there is 
a possibility to use the vacuum values of the corresponding LEC's to get an idea 
about the approximate order of their in-medium analogs. It seems that such 
a matching of the vacuum and in-medium LEC's could provide a quite impor- 
tant additional possibility allowing one to put some constraints on the values 
of the in-medium LEC's since the vast experimental data on the NN scattering 
in vacuum can be used in the calculations at finite (although moderate) densi- 
ties. It may also turn out useful when calculating the nuclear saturation curve 
since it helps to avoid the procedure of fixing the phenomenological in-medium 
amplitudes at different nuclear densities which would otherwise be needed to 
fix the density dependence of the effective couplings, making the whole EFT 
approach less predictable and reliable. Note that, although this matching can 
be carried out at any value of fi the choice /i = is rather inconvenient one 
since the values of the vacuum and in-medium LEC's compared at \i — are 
completely different so that the conclusion about the moderate density depen- 
dence of the effective couplings may no longer be true at /i = 0. Thus, the 
choice fi ~ ~ pp looks more natural. We emphasize that the observables, 
expressed via G-matrix do not depend on the value of fi chosen. However, it is 
much more difficult to relate the in-medium and vacuum NN-forces when the 
value fi = is taken for the subtraction point. 

We observe approximately 35% change in the value of Co(/i = m^^pp) 
when going from LO to NLO. It indicates that the chiral expansion is sys- 
tematic in a sense that the NLO corrections lead to the "NLO changes" of 
the effective couplings already determined at LO. The natural size of the in- 
medium effective range parameters, moderate changes experienced by the LO 
coupling constant Co and smallness of the NLO LEC's might, in principle, 
indicate the possibility of the perturbative calculations. However, in spite of 
this, it is still more useful to treat this problem in the nonperturbative man- 
ner. There are few reasons for this. Firstly, the overall (although distant) 
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goal of the EFT description is to derive both nuclear matter and the vacuum 
NN amplitude from the same Lagrangian. However, it is hard to say at what 
densities the dynamics becomes intrinsically nonperturbative, so it is better 
to treat the problem nonperturbatively from the beginning. The nonpertur- 
bative treatment may also turn out important to get the correct saturation 
curve since at some density lower than the normal nuclear one the scattering 
length starts departing from its natural value and some sort of the nonper- 
turbative approach becomes inevitable. Besides, the logarithmic terms in the 
chiral G-matrix are important to ensure correct threshold behavior of the in- 
medium NN-amplitude. Secondly, the NLO corrections themselves are quite 
significant. Thirdly, the nonperturbative treatment is required to study the 
phenomena of nuclear superfluidity p3 when the perturbation theory is not 
valid. Finally, in the processes involving both the nonzero density and tem- 
perature, such as heavy ion collisions, the value of the Fermi-momentum can 
effectively be lowered again making the nonperturbative treatment preferable. 

As we mentioned in the introduction , one of the most important (and yet 
unsolved) problems of constructing EFT at finite density is the formulation 
of the corresponding counting rules. The complete solution of this problem is 
possible only if pion effects and many-body forces are taken into account in a 
chirally invariant manner. However, as a first step in this direction one could 
formulate the "naive" counting rules for the LEC's Cqm) used in the effective 
Lagrangian. Since the main reason of the anomalous KSW counting C*2n ~ 
ma>"+ 1 su gg es ted by Kaplan et al. || is no longer the case in nuclear medium 
it would be natural to count like C 2n ~ ma^"+ 1 • However, one could modify 
this counting taking into account a further small scale Q ~ pp ~ ~ 
involved at nonzero density so that taking pp — > limit gives rise to the KSW 
scheme. So we count 



4-7T 

„ Till 



MA> n+1 



1 n+l 



-fx- (p F /n)(A-n) 



(17) 



Assuming fi ~ pp one indeed observes the moderate dependence of the LEC's 
on density. The above stated values of LEC's are consistent with this power 
counting. Note that in order to be fully consistent the density dependent nu- 
cleon mass should, in principle, be used in Eq.(17). However, to get the crude 
estimate of the effective couplings the use of vacuum value of nucleon mass 
is sufficient. It is clear that the exact form of the density dependence of the 
effective couplings is not unique. Any functional form providing the smooth 
interpolation between the "natural" counting C 2n ~ Mj\2n+ 1 &t normal nuclear 
density and the KSW counting C 2n ~ MA y +1 at zero density is acceptable. 
It is plausible however that all forms will give the same order-of-magnitude es- 
timate of the effective coupling. When calculating some observables at normal 
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nuclear density the simplest thing is to make use of the "natural counting". 
The further understanding of the density dependence of the in-medium LEC's 
could be achieved considering the nuclear saturation since the saturation curve 
would be quite sensitive to the exact values of the LEC's at different densi- 
ties. However, to treat the saturation in a consistent manner the higher partial 
waves (at least P and D waves) and pion effects must be included. We rel- 
egate the discussion of nuclear saturation to the separate paper. Note that 
although the value of the in-medium cutoff parameter A may also differ from 
its vacuum value we do not expect this difference to be significant. It seems to 
be rather unlikely that the introduction of the additional low-momentum scale 
Pf can lead to the large changes of the cutoff A, which is meant to be the in- 
direct manifestation of the truly short-range physics which is only moderately 
touched upon by the medium effects. 

Assuming that the values of the effective couplings are fixed at zero nu- 
clear density and using the above written expression for the couplings as the 
functions of density one can determine the LEC's at any intermediate value of 
the Fermi-momentum. 

Now the remark, concerning the naturalness criteria [I5| is in order. The 



naturalness concept in the context of nuclear interactions was elaborated in |7[] . 
According to this concept of an individual term in the effective Lagrangian can 
schematically be written as some dimensionless factor of order unity multiplied 
by the certain combination of scale factors, characterizing the scales involved 





1 


~d' 


[ f 2 A . 




A 



(M) 2 - (if 



Applying the scaling rules developed in [7[] to extract all scale factors and 
assuming that the cutoff parameter A ~ 500 MeV one finds 00(02) — 0(1). 
Thus the dimensionless coefficients are indeed compatible with naturalness. 
4. Low-density limit. Perturbative case 

The low density limit has often served as a somewhat toy many-body prob- 
lem, useful for understanding the structure of the theory and checking the 
consistency of the approximations made. In this section we consider the case 
when all interactions are natural so that the perturbative treatment is valid. 
The corresponding expression for the potential energy per particle of the dilute 
Fermi-gas is well known and was derived by traditional many-body methods 
long time ago. We show in this section that in the limit when the perturbative 
expansion is valid (G ~ V) the LO term of the low-density expansion can be 
reconstructed within our approach. Note that the low-density expansion in- 
cluding both leading and higher order terms has recently been rederived in ref. 
H in the framework of perturbative EFT. The pairing phenomena requiring 
the nonperturbative treatment will be considered in the next section. 
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The potential energy of the system of the interacting fermions is 



1 W 



(19) 



The summation goes over the states with momenta below pp. Since by as- 
sumption the perturbative expansion is valid we put (G ~ V ~ C ). The 
explicit expression for the S'-wave part of ground-state energy per particle is 



E 



3p2p + -^—(2T + 1) ( PF dkk 2 G{k,p F ) 

Jo 



A 10M ' vrV 

PF-k 



(PF-P 2 ) 1/2 
PF-P 



k 2 



2Pk 



(20) 



where P is the total momentum of the pair. Since we consider the low-density 
limit it is legitimate to use the relation Co = ^ between the effective coupling 
Co and scattering length a. After the momentum integration is done we obtain 
the following lowest order expression for the ratio 4 in the case of the dilute 
Fermi-system 



E 
A 



Mo 



3 

To 



l 

3^ 



(p F a) + 0{p 2 F ) 



(21) 



which coincides with the well known textbook expression ]T6] . The next ana- 
lytic in pf terms of the low-density expansion can be derived if NLO term with 
the effective coupling C2 is included. The result also agrees with that obtained 
by standard many-body technique. One notes that it may be nontrivial to 
go further in the low-density expansion since some higher order terms appear 
in the conventional expression multiplied by powers of ln{ppa). As argued in 
Ref. in order to reproduce these terms in the framework of the perturba- 
tive EFT the graphs describing three-to-three scattering amplitude must be 
included. The point is that the diagram for the two-to-two scattering contain 
only power divergences and thus, being regularized by means of dimensional 
regularization with minimal subtraction, can only lead to the terms which are 
analytic in pf/K.. Note that it may not be so in the nonperturbative situa- 
tion with unnatural scattering length since the G-matrix from Eq.(12) already 
contains the logarithms of pf needed to reproduce the terms with ln(pFd) in 
the low-density expansion so that no additional 3-body terms are needed. Of 
course, in the case of large scattering length the low-density expansion has 
a little practical use since it is valid for very low values of pp, satisfying the 
condition ppa < 1 

5. Low-density limit. Nonperturbative case 
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It is well known that nuclear matter is superfluid. The nucleon-nucleon 
interactions is attractive in the 1 Sq channel so that the corresponding ampli- 
tude develops the singularity, no matter how weak the interaction is. This 
singularity is cured by the formation of Cooper pairs. The superfluid prop- 
erties of nuclear matter are important in the study of neutron stars [O and 



heavy nuclei close to the drip line [18|. The studies of nuclear superfluidity 



have basically been carried out using pairing matrix element given by the bare 
nucleon-nucleon interactions as the elementary input entering the many-body 
calculations so it would be useful to see how nuclear superfluidity is described 
when the basic nucleon-nucleon interaction is treated in the framework of EFT. 

The gap exhibited by the system is formed at the Fermi surface. Since 
Cooper pair condensation cannot be described by means of a perturbation 
theory a nonperturbative treatment is required. The standard approaches to 
the problem of nuclear superfluidity are based on the traditional many-body 
methods such as coupled-cluster theory |1| or self-consistent Green function 



theory [£(]] . To our knowledge the first attempts to analyze the nuclear super- 
fluidity in the framework of EFT was made in Refs. f21||2~2|j . The discussion 
below follows the line which is somewhat similar to that in Ref . [^TJ . The aim 
of this section is to use the RG equations to estimate the size of the superfluid 
gap at small nuclear densities. In the RG technique the size of the gap is 
determined by the position of the singularity in a running effective coupling 



near the Fermi surface 23 . Let's define 



—arctan(-) = t , —arctani—)) = t. (22) 
Pf Pf Pf Pf 

After substitution to Eq.(16) one gets 

r (t n ) - C {t ,p F )27T 2 

Co ^ PF) -2^ + M PF (t -t)c (t , PF y (23) 

We take the limit p F — > so that the effective coupling Cq extracted from the 
nucleon-nucleon scattering in vacuum can be used 

47T 1 

We denote a the vacuum scattering length. One can see that the coupling Co 
has singularity at 

"* ** A. (25) 



M C (t) 7T 

The pairing gap is expected to be of size ~ exp^^u*). Assuming t = one 
can rewrite the expression for the gap in terms of the scattering length 
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A ~ bexp(-^ ). (26) 
2p F a 

This expression coincides with one obtained in Ref. [22]. However, as pointed 
out in || the choice \i — is not the optimal one. The choice fi ~ p makes the 
power counting more transparent. Moreover, at p ~ the density dependence 
of the Co may become nonnegligible even at low densities so the assumption 
Co(/x,Pf) — Co({i,pf = 0) is better justified when /x > pp. The scattering am- 
plitude is independent of the subtraction point so the choice of this point is just 
a matter of practical convenience. We used the value Co(/i) = — 2.23/m 2 de- 
termined at /i = X.lfmT 1 and calculated the gap A(p F ) as a function of Fermi 
momentum. One notes that the RG analysis gives the order of magnitude 
estimate of the gap but cannot provide the determination of the preexponen- 
tial factor which may be important at small pp. This factor can be found by 



matching to the known results, obtained at low density [16]. The result of this 
matching is 



PF 



(27) 



2M 

The results of calculations of A(p F ) are 

A(p F = O.l/m- 1 ) = 0.075MeV, A(p F = 0.2/m" 1 ) = 0.3MeV, 

A(p F = OAfm- 1 ) = 0.7QMeV, A(p F = 0.6/m" 1 ) = 2.39MeV, (28) 

We stopped calculations of the gap at p F = 0.6/m -1 since at higher values 
of Fermi momentum the condition /i > p F is no longer valid and the density 
dependence of the effective couplings may affect the results of calculations. 

Our results are in a qualitative agreement with the calculations using more 
traditional approaches. However, as well known from the phenomenological 
studies at p F ~ 1.7/m _1 the gap vanishes. The physical reason is that at 
such densities the attraction in the 1 Sq channel is replaced by the repulsion 
so that the formation of the fermion condensate is no longer possible. The 
EFT treatment must be modified to incorporate the effect of the decrease 
of a superfluid gap with increasing density since the ration p F /A is not a 
reliable expansion parameter any longer and density dependence of the LEC's 
should be taken into account. The problem of the EFT description of nuclear 
superfluid gap in the wider range of densities will be addressed in a separate 
paper. 

6. Nuclear matter observables 

Let's now discuss the results of calculations of potential energy of nuclear 
matter. As in the low-density case the Eq.(20) is used. Let's start from 
the results of calculations at normal nuclear density po = 1.36/m -1 . In the 
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LO perturbative calculation assuming G ~ C one gets ~ -ll.VMeV. 

Making use of the nonperturbative expression for the G - matrix but keeping 
only LO in the chiral expansion of the effective nucleon-nucleon potential gives 
rise to the value ~ —17.0MeV. The inclusion of the NLO terms in 

the chiral Lagrangian leads to ^ s °) ~ —13.0MeV . So the corresponding 
correction is about 30% and can be attributed to NLO. One notes that both 
^gg) an( j lo effective coupling C experience NLO corrections when NLO 
terms are included in the effective Lagrangian. It justifies one more time the 
statement about the consistency of the corresponding chiral expansion. The 
similar calculations performed in the triplet S-wave channel leads to the value 
~ -17.3(-13.2)MeV in LO (NLO). However, there is an important 
difference between the singlet and triplet channels because of the S-D mixing 
arising at NLO in the triplet channel. This effect is not considered in this 
paper since it requires the inclusion of the explicit pion degrees of freedom 
in the Lagrangian. The issue of the pionic effects is discussed in more details 
below. Here we only note that the pion effects must eventually be included in a 
chirally symmetric manner consistent with the in-medium counting rules that 
are yet to be established. It makes the whole issue of the pionic effects quite 
complicated. To calculate the value of ^ at lower densities we use Eq.(17) to 
fix the effective couplings Co and C 2 - We use the value A = 2.05/m~ 1 for 
the effective cutoff parameter. The effective couplings determined at different 
densities and at ji = m w are 

C (fx,p F = 1.2/m" 1 ) = -2Alfm 2 ;C 2 (fi,p F = 1.2/m" 1 ) = 0.98/m 4 
C ((i,Pf = 1.0/m -1 ) = -2A8fm 2 ;C 2 (fi,PF = 1.0/m -1 ) = 1.05/m 4 
C (h,Pf = 0.8 fm- 1 ) = -2.53/m 2 ; C 2 (fx,p F = O^fm' 1 ) = 1.16/m 4 (29) 

For the potential energy per particle we find ^ = - 11.5 MeV (-9.7 MeV, - 8.1 
MeV) at pf = 12 /to -1 (1.0/m -1 , 0.8/m -1 ) respectively. These results look 
quite reasonable, although somewhat smaller than the values obtained in more 
traditional approaches. We stopped the calculations at pp = 1.36/m _1 which 
is rather close to the nuclear saturation point. The approach should be mod- 
ified in the several aspects in order to be applied to the nuclear systems with 
the density higher then the normal nuclear one. First, the Fermi-momentum, 
being only "marginally light scale" already at normal nuclear density becomes 
closer to the generic "heavy" scale. Consequently, the counting rules must 
be modified and the value of the effective cutoff increased. At some critical 
density the value of cutoff can reach the point, where the other degrees of 
freedom such as vector mesons and A-isobars are to be included explicitly. It 
brings in the other sources of uncertainties which are rather difficult to control. 
Moreover, whereas at normal nuclear density and typical energy scale Q ~ 100 
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MeV the relativistic corrections enter at higher orders, at some critical den- 
sity the whole formalism should be formulated in a fully relativistic way. It 
considerably complicates the practical use of EFT. The counting rules must 
be significantly modified to take into account the excitations both the Fermi 
and Dirac seas and the fact that the nucleon mass can no longer be considered 
as a genuine "heavy" scale. On the other hand, at densities much higher than 
the normal nuclear one the matter was shown to be a color superconductor 



24] where quark and gluon degrees of freedom must be considered explicitly. 



When the density is lowered a color superconducting phase should match the 
EFT, formulated in terms of baryons and mesons. Carring out this matching 
is a highly nontrivial problem, requiring at least qualitative understanding of 
the mechanism of hadronization at high densities |[25|| . Therefore some model 
of how confinement occurs at high densities is probably needed. 

As already mentioned above, to be completely consistent, many-body forces 
and pion effects must be included in the EFT treatment. The contribution from 
the many-body forces is expected to be rather small. There are at least two 
reasons for such an expectations. First, each additional interacting nucleon 
line will bring in the additional small factor Q 3 ~ p 3 F in the expression for 
the binding energy of nuclear matter. Second, following the line of arguments 
suggested in [^6| the contributions of the N-nucleon forces can schematically 
be written as Vn ~ Q N /A^ -1 where A is the generic large-mass QCD scale 
close to 1 GeV and Q is the typical low-energy scale taken to be about 0.1 
GeV. Using these numbers one can see that the contribution of the three-body 
forces is order of magnitude smaller compared to the two body ones and that 
four-body interaction can safely be neglected. We stress, however, that taking 
into account the 3-nucleon forces may turn out important to get the saturation 
point close to the empirical one. The expected smallness of the many-body 
forces compared to the two-body ones is also confirmed by the results of several 
other studies both for nuclear matter and for finite nuclei with A > 3 |p7|| . 

Pion degrees of freedom, being another important ingredient of any nuclear 
EFT, seem to average to NLO effects, when the S-wave part of nuclear matter 
binding energy is being calculated. One notes that this conclusion is partly 
supported by the results, obtained in the phenomenologically very successful 
Skyrme approach [ |2"5| , |2T| . The Skyrme interaction has the pointlike form and 
does not contain the explicit pions. In some sense the Skyrme interaction can 
be considered as a parameterization of the G- matrix [30|, written in the form 
closed to that we used. Furthermore, as shown in |3l[] adding explicit pions 
can lead only to the moderate changes of the effective couplings. This fact can 
qualitatively be explained by the partial cancelation from the pointlike inter- 
action and one-pion-exchange interaction. The existence of such a cancelation 
was also pointed out in H. However, it is worth noting that the relative con- 
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tribution due the pionic effects may depend on the partial channel considered. 
For example, we expect that the pion degrees of freedom are more important 
in the P-wave part of nuclear potential energy that in the S-wave one. 

There are a few other important issues worth discussing. First, the cal- 
culations must be self-consistent in a sense that both G-matrix and effective 
nucleon potential, entering the BG equation should be determined at the same 
order from the system of coupled equations since effective nucleon potential 
is itself determined by the G-matrix. This is not done yet. Second remark 
concerns the problem of removing of the off-shell ambiguities. EFT offers, in 
principle, the potential possibility to avoid this difficulty order by order. In 
practice, however, it may turn out that the cancelation of the unwanted contri- 
butions can be achieved only if all graphs, relevant at a given chiral order are 
included. It may turn out very complicated technical task. Third, the nuclear 
dressing of pions is to be considered. In the other words, pion self-energy must 
be computed to the order dictated by the counting rules. 

In conclusion let us list the main ingredients which have to be included in 
any consistent nuclear EFT. 

1. The local interaction terms as well as pion degrees of freedom must be 
included, probably up to NLO in the S, P and D-waves. 

2. Nucleon and pion self-energies should be computed at a given chiral 
order. 

3. To explain the saturation the three- nucleon forces seem to be needed. 

4. The relations between the in-medium effective couplings and those aris- 
ing in the vacuum 2-body and few-body calculations have to be established. 

5. The reasonable way to generalize the nuclear matter EFT to treat the 
finite nuclei should be found. 

6. All that must be supplied with the consistent formulation of the in- 
medium counting rules. This is probably a crucial issue for future studies. 

As seen from this list we are at the beginning of the road. Many things 
must be done to formulate the entirely consistent and yet practical approach. 
However, the results already obtained on this way give a good reason to hope 
on a further progress in applying of the EFT methods to nuclear physics. 



16 



REFERENCES 



[1] A. Manohar, Effective Field Theories, |hep-ph/9506222| ; D. B. Kaplan, Effec- 



tive Field Theories, nucl-th/9506035 



S. Weinberg, Nucl. Phys. B363, (1991) 3. 

D. B. Kaplan, M. Savage, and M. B. Wise, Nucl. Phys., B534, (1998) 329 
M. C. Birse, J. A. McGovern and K. G. Richardson, Phys. Lett., B464, (1999) 
169. 



B. V. Krippa, |nucl-th/0002019 



H.-W. Hammer and R. Furnstahl, |nucl-th/0004043 . 

R. Furnstahl, B. Serot and H.-B. Tang, Nucl. Phys., A618, (1997) 445. 
M. Lutz, Nucl. Phys., A642, (1998) 171. 

A. Bulgac, G . A. Miller and M. Strikman, Phys. Rev. C56, (1997) 3307. 
M.Jaminon and C. Mahaux, Phys.Rev. C40, (1989) 354. 

S. Beane, T. D. Cohen and D. Phillips, Nucl. Phys., A632, (1997) 445. 

B. V. Krippa and B. L. G. Bakker, |hep-ph/0004178| ; 
J, Gegelia Phys. Lett., B429, (1998) 227. 

A. B. Migdal, Sov. Phys. JETP 10 (1960) 176. 

H. Georgi and A. Manohar, Nucl. Phys., B234, (1984) 189. 

A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems 

(McGraw-Hill, New- York, 1991) 

C. J. Pethick and D. G. Ravenhall, Annu. Rev. Nucl. Part. Sci. 45 (1995) 429. 
A. C. Muller and B. M. Sherill, Annu. Rev. Nucl. Part. Sci. 43 (1993) 529. 
R. F. Bishop, Theor. Chim. Acta 80 (1991) 95. 

W. H. Dickhoff and H. Muther, Rep. Prog. Phys. 55 (1992) 1947 

J. Hormuzdiar and S. Hsu, nucl-th/9811017 

T. Papenbrock and G. F. Bertsch, Phys.Rev. C59, (1999) 2052. 

J. Polchinski, in Proc. of the 1992 TASI, eds J. Harvey and J. Polchinski (World 

Scientific, Singapore 1993). 

M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett B422, (1998) 247; R.Rapp, 
T. Schafer, E. V. Shyryak and M. Velkovsky, Phys. Rev. Lett. 81, (1998) 53. 
S.Pepin, M. C. Birse and J. McGovern, |hep-ph/9912475 



J. Friar, Few-Body Systems 22, (1997) 161 

J.Carlson and R. Chiavilla, Rev. Mod. Phys. 70 (1998) 743 

T. H. R. Skyrme, Nucl. Phys., 9, (1959) 615. 

R. Furnstahl and J. C. Hackworth, Phys.Rev. C56 (1997) 2875. 

D. Vautherin and D. Brink, Phys.Rev. C5, (1972) 629. 

B. Lynn, Nucl. Phys, B402, (1993) 281 



17 



